Enhancement of the critical temperature in iron-pnictide superconductors by finite 

size effects 
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Recent experiments have shown that, in agreement with previous theoretical predictions, super- 
conductivity in metallic nanostructures can be enhanced with respect to the bulk (L — > oo) limit. 
Motivated by these results we study finite size effects (FSE) in an iron-pnictide superconductor. For 
realistic values of the bulk critical temperature jj ulk ~ 20 — 50K, we find that, in the nanoscale 
region L ~ 10 nm, T C (L) has a complicated oscillating pattern as a function of the system size 
L. A substantial enhancement of T c with respect to the bulk limit is observed for different bound- 
ary conditions, geometries and two microscopic models of superconductivity. Thermal fluctuations, 
which break long range order, are still small in this region. Finally we show that the differential 
conductance, an experimental observable, is also very sensitive to FSE. 

PACS numbers: 74.70.Xa,74.78.Na,74.20.Fg,74.78.-w 
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Recent technological developments 0, Q have opened 
the possibility for the synthesis and characterization of 
high quality nano-structures, a recurrent problem in pre- 
vious [3j experimental studies of superconductivity in the 
nanoscale region. As an example, it is now possible Q 
to track experimentally the evolution of superconductiv- 
ity in single, clean and isolated metallic nanoparticles of 
typical size L > 5 nm by using low temperature scanning 
tunnelling microscopy /spectroscopy. 
These advances have revived the interest in nanoscale su- 
perconductors. Theoretical studies on this problem were 
pioneered by Anderson [j| more than fifty years ago. 
At zero temperature the exact Richardson's equations 
@ provides a satisfactory account of quantum fluctua- 
tions and other deviations from mean-field behaviour in 
nanosuperconductors. At finite temperature the role of 
thermal fluctuations is well described by the path inte- 
gral formalism Q. In the limit of negligible disorder, 
FSE in conventional superconductors have been studied 
for different systems: an harmonic oscillator [8|]._rect an- 
gular nanofilms 
a chaotic grain 



a nanowire [9( , a sphere 12] , and 
13l Il4 | . The main conclusions of these 



studies are that: a) FSE can enhance or suppress T c sub- 
stantially provided that the superconducting coherence 
length £ is not much smaller than L, b) FSE are stronger 
the smaller and more symmetric the grain is, c) thermal 
fluctuations are controlled by the parameter 7 = f3y/S/T c 
Q where S is the mean level spacing around the Fermi 
energy and (3 » 1/2. For 7 < 1 long-range order still 
holds as only a narrow region ~ "/T c around T c is af- 
fected by thermal fluctuations. These conditions set the 
minimum L for which enhancement of T c is observed. 
A natural question to ask, especially after the recent dis- 
covery of iron pnictides |15[ , is whether these results are 
of relevance for high T c superconductors. To the best of 
our knowledge not much is known about FSE in high T c 
materials. Cuprate nanowires [l6| and carbon nanotubes 
17] have been recently investigated experimentally for 



sizes where thermal and quantum fluctuations are dom- 
inant. We are not aware of theoretical studies of FSE 
in high T c superconductors (for a recent study on multi- 
band superconductors see [la]). This paper is a first step 
in this direction. 

We investigate superconductivity in clean, rectangular 
nanofilms of iron pnictides. All our calculations are in 
two dimensional geometries as it is believed that this is 
the effective dimensionality relevant for superconductiv- 
ity in these materials. Iron-pnictides are more suitable 
than cuprates because, unlike cuprates, superconductiv- 
ity in iron pnictides seems to be well described by a mean 
field approach based on Fermi liquid theory. This is a 
clear advantage as the ideas and techniques introduced 
in conventional nanosuperconductors are still applicable. 
Moreover FSE are easier to detect in iron-pnictides as its 
£ is larger than that of cuprates. 

Models: We study FSE in two popular microscopic 
models of iron based superconductors: the minimal two 
band model of [l9[ and the five band model of [2(j. It 
is important to investigate FSE in more than one model 
as, in these materials, there is not yet consensus about 
fundamental issues such as the mechanism for super- 
conductivity (recent experiments [2l| suggest spin ex- 
change interactions) or the leading pairing symmetry (see 
[H El E H| for different proposals). 
First we introduce the two band model of [l9j: the dis- 
persion relation in the hole and electron band is given 



by, e h (k) = t^(cosk x 



e e (k 
ing 



cos ky) +ttj; cos k x cos k y 



cos k x + cos k y ) + 1\ cos % cos 



e and 
e e . Follow- 



22J we choose (ti,i 2 ,e) to be (0.30,0.24,-0.6) 
for the hole band and (1.14,0.74,1.70) for the electron 
band. The bulk density of states is N h {0) = 0.74/eV 
and N e (0) = 0.285/eV. Our results are robust to small 
changes of these parameters. Superconductivity can be 
investigated by a mean field approach with two order pa- 
rameters, Ah and A e , related by the following self con- 
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FIG. 1. <5T C (L)/T c bulk (typical deviation of T C (L) from the 
bulk limit T c bulk = T C {L — > oo) in a square film of side L 
for different bulk T c and PBC. Upper: two band model fTJ. 
Lower: five band model [20|. Results of both models are qual- 
itatively similar. Typical enhancement (suppression) of T c of 
more than 50% is observed for L ~ 10 nm. 
Inset: Critical temperature T C (L) as a function of L for 
T bulk « 30K. 

sistent conditions: 

tanh(Ml) tanh(^)) 
A ' 1 - ~ y VhhAh 2E h (k') + VheAe 2E e (k>) ' 

tanh(^l) tanh(^l) 
Ae - ~ y VehAh 2E h (k>) + VeeAe 2E e {k') 

(1) 

where, Vh e ,he,V eei hh are the inter-band and intra- 
band coupling constants respectively and Eh, e (k) — 

^J(^h,e(k) — + A 2 e (k) are the energies of the quasi- 
particle excitations and fi is the chemical potential. The 
sums above are restricted to energies ljafm = 100 meV 
around fi. FSE are studied by simply solving (JT|) for a 
rectangular geometry. We note that the different cou- 
pling constants are in principle momentum dependent 
Vi = Vi(k,k'). In [19] it was indeed employed a mo- 
mentum dependent spin exchange interaction. It was 
however found that the order parameters were weakly 
dependent on momentum. Based on this result we have 
neglected any momentum dependence of the bulk inter- 
action. In finite size systems the coupling constant gets 
a momentum/energy dependence [9l. lloL Il3| even if the 



bulk interaction is constant. Assuming a contact in- 
teraction, it is easy to show 0, [3 that V, oc I(k, k') 
where I(k, k') are the so called matrix elements I(k, k') cx 
L 2 J dxdyip^(x, y)i\)\i (x, y) and ipk are the eigenvectors of 
the one-body problem. In the case of chaotic grains [13| , 
or rectangular grains with Dirichlet boundary conditions 
[Toj . this correction induces a non-oscillatory enhance- 
ment of the energy gap in the nanoscale region (see inset 
Fig. [2] (lower)). More specifically, the coupling constant 
effectively increases as L decreases 0, [T(Jll3j]. It is also 
clear from the definition of I(k, k') that, in a rectangular 
system with periodic boundary conditions (PBC), the in- 
tegral above is momentum independent. Due to this ad- 
ditional simplification we have decided to use PBC. We 
note that this is also a sensible choice to describe granu- 
lar systems, where many grains are strongly coupled, or 
superconducting rings. Moreover results for other bound- 
ary conditions are qualitatively similar (see Fig. [2]). 
We also investigate the five band model of which pro- 
vides a fully quantitative description of the band struc- 
ture of an iron based superconductor within 2 eV of the 
Fermi energy. For details of the model we refer to the 
appendix of [2(]j ]. Superconductivity, assumed to be of 
a s ± type A(k) = A(L) cos(k x ) cos(k y ) , is studied in a 
BCS mean field fashion involving the five bands. FSE arc 
investigated by solving exactly the mean field gap equa- 
tions for a given geometry and loafm — 80 meV. 
Results: Our first task is to compute T C (L) (see inset 
Fig. QJ numerically for the two models above. Calcula- 
tions are restricted to the range L £ [10, 40]nm in order to 
keep the effect of thermal fluctuations small. For each L 
the chemical potential [i was computed keeping constant 
the electron density. For the sake of simplicity we present 
results for a square film with PBC. It will be shown that 
qualitatively similar results are also obtained for rectan- 
gular films and other boundary conditions (see Fig. [2]). 
The value of the coupling constants are set to lead to 
typical values of T^ ulk £ [30, 50]K in iron based super- 
conductors. FSE are rather insensitive to small changes 
of these values provided that T^ ulk is still the same (see 
@). For a more quantitative assessment of the typical 
maximum enhancement of T c we also plot 8T C (L") /T~ , 
the typical deviation of T C (L) in units of T c bulk . ST C (L) is 
obtained by evaluating the standard deviation of T C (L) 
in consecutive ~ 3 nm intervals where T C {L) is computed 
in steps of 0.04nm. In Fig. [1] we observe that: a) the 
typical deviations with respect to the T c bulk is quite large 
(> 50%) in the L ~ 10 nm region, b) even for L ps 10 nm 
and T c ps 30K thermal fluctuations are still small since 
7 ~ 0.3 < 1. We note that in our case 5 is the geometrical 
mean (see ^) of the mean level spacing in each band, b) 
T C (L) is an oscillating function of L with local maxima 
T C (L) much higher than in the Xjuik ij m ^ ( see inset), c) 
FSE increase as L or T c bulk decreases, d) results for both 
models are similar. 

In order to understand the origin of these features we 
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FIG. 2. Upper: <5T C (L)/T c bulk from JTJ) for a square of side 
L and two rectangular films (L x = aL,L y = L/a) with a — 
1.618034 (rectangle 1), a — 1.2334 (rectangle 2) respectively 
and Tc Ulk ~ 32 K. The intricate band structure and pairing 
pattern make the dependence on the shape weaker than in 
single band superconductors with parabolic spectrum. Inset: 
T C (L) for the square and rectangle. Lower: ST C (L) /T bulk in a 
square film, T c bulk « 32K, for Dirichlet (including the effect of 
the matrix elements I(k, k')), and PBC. In the latter, relative 
deviations are stronger due to the additional level degeneracy 
but absolute deviations for small sizes can be still stronger 
for Dirichlet due to the effect of I(k,k') (see text). Inset: 
Comparison between T C (L) for Dirichlet and PBC. The solid 
lines stand for the local T C (L) average in the same intervals 
(~ 3nm) in which 5T C (£)/T c bulk is evaluated. 



study ([IJ in the limit V ee — Vhh = (this is a reason- 
able assumption as inter-band pair hopping is known to 
be the dominant mechanism) where an analytical estima- 
tion of T C (L) is feasible by using semiclassical techniques 
I13L 1 24| . For single band superconductors it was shown 
in [13J that, in the region in which mean-held techniques 
are applicable, non-monotonic deviations from the bulk 
limit are well described by simply replacing sums by inte- 
grals in the gap equation, ([lj in our case, and extracting 
from the integral a modified DOS A^(0). Technically this 
modified DOS is expressed as a sum over classical peri- 
odic orbits of the grain with lengths smaller than £. An 
explicit expression for a rectangular film can be found 
in [24J and in the Appendix of the second reference of 
[l3j . In practical terms a similar result is obtained by 
simply smoothing the DOS N(0) ~ J2 i 5(e — ej) over an 
energy scale hvp/n^ where vp is the Fermi velocity. It 
is straightforward to show that a similar approach in our 



case leads to, 

T C (L) « 1.136WAFM exp(-l/A cff ) (2) 



where A 



off 



= A 



bulk 



jV| e (O)jV^ (0) 



with Abu 



lk — 



v / V eh N e (0)V he N h (0), and £ fc(e ) stands for £ in the hole 
(electron) band. 

Several comments are in order: a) fluctuations in T C (L) 
are caused by finite size, not thermal, fluctuations of 
the spectral density around the Fermi energy. In sys- 
tems with a parabolic spectrum it is well known that 
level degeneracy induced by geometrical symmetries [24 1 
of the film (cube, sphere...) increase these fluctuations, 
b) (small) changes in A e ff, due to FSE, induce (exponen- 
tially) large changes in T c . Therefore small changes in 
the number of levels around the Fermi energy can lead 
to large changes in T c , c) FSE are strongly suppressed for 
(/L « 1 since, in this limit, the discrete nature of the 
spectrum is completely smoothed out and A e ff ~ Abuik- 
Since £ decreases as Abuik increases, deviations from the 
bulk limit for a given size L decrease as T^ ulk increases. 
This is what is observed in Fig. [1] 

The rest of predictions are now tested numerically for the 
two-band model ([1]). More specifically we study the de- 
pendence of FSE on boundary conditions and the shape 
of the grain. In Fig. [5] (upper) we compare the typ- 
ical deviation 8T C (L)/T^ ulk for two different rectangles 
and a square of the same total area. Deviations from 
the bulk limit are slightly larger in the square. This 
is due to the existence of a symmetry k x — > k y in the 
spectrum of the square which is not present in the rect- 
angle. As was mentioned previously, level degeneracy in 
the spectrum leads to stronger fluctuation in the number 
of levels around the Fermi energy and, according to ([2]) , 
to larger deviations of T C (L) from T^ ulk . We note how- 
ever that, the difference between the two geometries is 
much smaller than in single band superconductors with 
parabolic spectrum. The reason for that is that the level 
degeneracy of a parabolic spectrum is much larger than 
the one (k x — >• k y symmetry) present in the models we 
study. In Fig. d] (lower) we compare <5T C (L)/T c bulk for 
a square with PBC and Dirichlet boundary conditions 
including the effect of the matrix elements I(k,k'). It 
is observed that, for Dirichlet boundary conditions, rela- 
tive fluctuations are smaller. This is due to the additional 
level degeneracy for PBC where eigenstates with oppo- 
site momentum have the same eigenvalue. By relative 
we mean that, as in the other cases, the standard vari- 
ance is computed with respect to the local average T c , 
not the . As can be observed in the inset of Fig. [2] 
(upper), the overall effect of I(k,k') is to enhance signif- 
icantly [nUm the average T C {L). As a result, the abso- 
lute deviation from the bulk limit for Dirichlet boundary 
conditions is of the order or even larger than for PBC. 
In summary, in the range of T^ ulk £ [20, 50)K typical 
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FIG. 3. Differential conductance G(V) (|3} as a function of 
the bias voltage V for T « T c /3, T c » 31 K and different 
sizes. G(V) is rescaled by its value, Go(V), at T « 2T C . 
Left: Two band model |T}. The peak in G(V) coincides with 
mw(|A e |, |A ft |) in ©. Right: Five band model The 
peak closer to V = also corresponds with the value of the 
energy gap. The additional structure is a consequence of the 
combined effect of the band structure and level degeneracy. 
The substantial differences observed for small changes in L 
suggest that it is experimentally feasible to detect FSE by 
scanning tunnelling microscope techniques (STM). 



prediction for £ cannot be very different from the experi- 
mental one. Otherwise the experimental enhancement of 
FSE, which is controlled by the ratio £/L, will not occur 
in the range of sizes we are proposing. We have checked 
for the model of [2(| that the averaged £ over the Fermi 
surface is consistent with the experimental one. 

In summary, we have investigated FSE in supercon- 
ducting iron-pnictidcs nanostructures. Within a mean 
field approach we have identified a region L ~ 10 nm 
in which FSE can enhance T c substantially and thermal 
fluctuations, which break long range order, are not 
important. Qualitatively similar results are obtained 
for different boundary conditions, system geometries 
and band structure. Experimental observation of these 
effects is in principle feasible as observables such as the 
differential conductance are very sensitive to FSE. 
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of iron-pnictides materials, substantial deviations from 
the bulk limit are observed for different models, shapes 
and boundary conditions for sizes L ~ lOnm for which 
thermal fluctuations are not yet important. 

Finally we discuss to what extent is possible to observe 
experimentally FSE. In metallic nanograins the normal- 
ized differential conductance, 



G(V) = 



ATkr 



duj N s (u>) 



cosh 



-V 



2k B T 



(3) 



(where N s (cj) stands for the superconducting DOS) is 
one of the most popular experimental observables as it 
can be measured by STM techniques. In Fig. [3] we de- 
pict G{V) as a function of the bias voltage V in the two 
models above for different, but similar sizes. It is clearly 
observed that, in the nanoscale region, G(V) is highly 
sensitive to both FSE and the model used. Small changes 
in the system size induce strong modifications in G(V). 
This is another indication that it is in principle possi- 
ble to observe and study experimentally FSE by STM 
techniques. More precisely we believe that, in order to 
compare our results with experiments, the following con- 
ditions must be met: a) Fermi liquid theory and mean- 
field techniques must be applicable, b) it is technically 
feasible to manufacture single grains or granular materi- 
als of size L ~ 10 nm for which thermal fluctuations are 
negligible but FSE are still relevant, c) the theoretical 
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